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We study the interaction between Venus and the solar wind using a global three-dimensional 
self-consistent quasi-neutral hybrid (QNH) model. The model treats ions (H + , + ) as particles 
and electrons as a massless charge neutralising fluid. In the analysed Parker spiral interplanetary 
magnetic field (IMF) case (IMF = [8.09, 5.88, 0] nT) a notable north-south asymmetry of the 
magnetic field and plasma exists, especially in the properties of escaping planetary + ions. The 
asymmetry is associated with ion finite gyroradius effects. Furthermore, the IMF x-component 
results in a dawn-dusk asymmetry. Overall, the QNH model is found to reproduce the main observed 
plasma and magnetic field regions (the bow shock, the magnetosheath, the magnetic barrier and 
the magnetotail), implying the potential of the developed model to study the Venusian plasma 
environment and especially the non-thermal ion escape. 



I. INTRODUCTION 

Neither Venus nor Mars has a strong global intrinsic 
magnetic field and therefore the solar wind is able to 
flow close to the planets in regions where the exospheric 
neutral density is high. Because of the resulting direct 
interaction between the exosphere and the solar wind, 
ionized atmospheric neutrals are accelerated (picked up) 
by the solar wind electric field. Charge exchange reac- 
tions between the solar wind protons and planetary neu- 
trals also produce energetic neutral hydrogen atoms (H- 
ENAs). Observation of H-ENAs is thus a manifestation 
of the direct interaction between the solar wind plasma 
and planetary neutrals. Picked-up planetary + ions 
can also form energetic neutral oxygen atoms (O-ENA) 
via charge exchange process. The ion escape, H-ENAs, 
O-ENAs and electrons will be investigated at Venus and 
Mars by two identical instruments: ASPERA-4 on Venus 
Express [seeQ] and ASPERA-3 on Mars Express [see 01 . 

The Venus-solar wind interaction can be analysed self- 
consistently by three-dimensional (3-D) magnetohydro- 
dynamic (MHD) and quasi-neutral hybrid (QNH) mod- 
els. The MHD models [see, for example, 6, and references 
therein] are computationally cheaper than the QNH mod- 
els [see, for example. 13. fl3l| . thus allowing a bett er g rid 
resolution. Two-dimensional (2-D) QNH models [I3 al- 
low similar or even better resolution than 3-D MHD mod- 
els, but then e.g. the 3-D draping of the magnetic field 
around the planet cannot be efficiently studied. An ad- 
vantage of a 3-D multi-ion species QNH model, as com- 
pared with the 3-D multifluid MHD model, is that it 
includes finite ion gyroradius effects. Furthermore, dif- 
ferent ion species are allowed to have different velocities 
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and temperatures. Even arbitrary non-Maxwellian dis- 
tributions functions are allowed. Therefore, the QNH 
model makes it possible to increase our understanding of 
those plasma environments where the ion gyroradius is 
comparable to or larger than the planet size. 

In this paper the Venusian plasma environment is stud- 
ied with a recently developed 3-D QNH model. The 
model is based on previous 3-D QNH models that have 
been used to study plasma and magnetic fields near Mars 
Mercury yj, Saturnian moon Titan [9] and the Moon 
. The analysis focuses on the asymmetries caused 
by kinetic effects, the role of the IMF (interplanetary 
magnetic field) x-component and the properties of the 
escaping planetary + ions. 

The paper is organized as follows. First the developed 
Venus model is described. Then the basic properties of 
H + ions, + ions and the magnetic field near Venus are 
presented. 



II. MODEL DESCRIPTION 

A. Quasi-neutral hybrid model (QNH) 

The quasi-neutral hybrid model (QNH) is based on 
the QNH models that have been used earlier to study 
Mars, Mercury, Titan, and the Moon 0, H H ED] . Here 
we describe only some of its basic features and its new 
possibilities. 

In this work, the coordinate system is a Cartesian sys- 
tem with the origin fixed at the centre of Venus. The 
x-axis points from Venus towards the Sun. The z-axis 
is defined by u z = u x x v/\u x x v|, where u x is a unit 
vector in the positive x-direction, and v is the orbital 
velocity of Venus around the Sun, i.e., z is perpendicular 
to the orbital plane of Venus. The y-sxis completes the 
right-handed coordinate system. 

Figure G] illustrates the coordinate system, grid struc- 
ture and terminology employed. The size of the Simula- 
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Figure 1: Illustration of the used simulation box and the 
adopted terminology. See text for details. 



tion box is — 3Ry < x < 3Ry and — ARy < y,z < ARy 
and the radius of the obstacle to the plasma flow is taken 
to be Ry = 6051.8 km. The grid size is one tenth of 
the obstacle radius (605 km). The upstream plasma pa- 
rameters analysed in this paper are nsw = 14 cm -3 , 
Usw = [-430, 0, 0] km s _1 and the thermal velocity 50 
km s _1 . The IMF corresponds to a Parker spiral angle 
(= arctan(^B2 + B 2 J\B X \) of 36°: B sw = [cos(36°), 
sin(36°), 0] 10 nT = [8.09, 5.88, 0] nT. 

It should be noted that the aberration caused by the 
orbital motion of Venus of about 35 km s _1 caused by 
the rotation of the planet around the Sun had not been 
taken into account. In the used solar wind velocity of 430 
km s _1 the aberration angle would have been about 5° 
(= arctan(35 km s _1 / 430 km s -1 )). Had the aberration 
been taken into account, the solution presented in this 
paper would have been rotated about 5° clockwise around 
the +z-axis if the results were presented in a coordinate 
system where x-axis points from Venus towards the Sun. 
In the rotation, the IMF x and y components would also 
have been slightly changed. 

As depicted in Fig. GJ in this paper the z > hemi- 
sphere is referred to as the +Esw hemisphere. The con- 
vective electric field in the solar wind Esw (= —Usw x 
Bsw) on the XY plane points toward the +Esw hemi- 
sphere. The opposite hemisphere is referred to as the 
-Esw hemisphere. The terminology reminds us of the 
fact that the results presented depend on the the choice 



of the direction of the IMF and that the asymmetries ro- 
tate around the x-axis along with the IMF. The +Esw 
and -Esw hemispheres are important when the properties 
of the + ions are considered, as will be seen in the next 
section. 

In addition, the y < hemisphere is called the BS perp 
hemisphere because there are points on the bow shock 
on the y < hemisphere where the IMF is perpendicular 
and quasi-perpendicular to the bow shock in the anal- 
ysed IMF case. Similarly, the opposite y > hemisphere 
is referred to as the BS par because there are points on the 
bow shock on the y > hemisphere where the IMF is par- 
allel and quasi-parallel to the bow shock in the analysed 
IMF case. It should be noted that if one had analysed 
the IMF B x < case the BS per p hemisphere would have 
been at the y > hemisphere. The BS pe rp/BS par hemi- 
spheres will be discussed in the next section when the 
asymmetries in the magnetosheath and the magnetic tail 
lobes are considered. 

The model contains two ion species, H + and + . The 
ions are modelled as particles and they are accelerated 
by the Lorentz force: 

mH+ ^ = e ( E + v H+ xB ) (!) 
m ° + ^dF = e ( E + v o + xB )- ( 2 ) 

Here v, m, E, B and e are the ion velocity, ion mass, elec- 
tric field, magnetic field and elementary charge, respec- 
tively. The subscript H + (0 + ) refers to protons (oxygen 
ions). 

Electrons form a massless fluid and their equation of 
motion (the momentum equation) is 

E + U e xB = 7? a J, (3) 

where n e , U e , J and r] a are the electron density, electron 
bulk velocity, electric current density and the anoma- 
lous resistivity, respectively. In the hybrid model the 
resistivity function can be a fully 3-D function, i.e.,7? a = 
rj a (x,y,z). At Venus the resistivity is largest near the 
planet, where the electron-neutral collision frequency be- 
comes high, and smallest in the undisturbed solar wind. 
In this paper, however, no ad hoc 3-D or 2-D resistivity 
models are developed. In contrast, a constant resistivity 
230 kf^m is used everywhere outside Venus for simplicity. 
The electric current is defined as 

J = e(n H +V H + + n +V + - n e U e ). (4) 

Here rii and V$ are the number density and average veloc- 
ity of the ith ion species. The electric current density J 
is related to the magnetic field B according to Ampere's 
law: 

VxB = fioJ, (5) 
where /jLq is the vacuum permeability. 



The QNH model assumes quasi-neutrality, i.e., 

n e =n u ++n +. (6) 
The magnetic field is propagated by Faraday's law 
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where the electric field is derived from Eq. 03 

The solution of Equations G] - proceeds as follows. 
First, the total current density is derived from the mag- 
netic field using Ampere's law ©. Second, the velocity 
field of the electron fluid is derived from the electric cur- 
rent, the number density and the average speed of ions 
by Equations El and El Third, the updated electric field 
is derived from the updated magnetic field, the current 
density and the velocity field of the electron fluid by us- 
ing Equation Finally, the ions are moved according to 
Equations G] and [21 by using the new electric and magnetic 
fields. The movement of ions results in updated plasma 
densities and bulk velocities and the magnetic field can 
be propagated by using the updated electric field in Fara- 
day's law (0) [seeH, for the details of the used algorithms]. 

The obstacle within the simulation box is assumed to 
represent an ideally conducting ball inside which the re- 
sistivity is zero. We apply fully absorbing boundary con- 
ditions on the surface of the obstacle by taking an ion 
away from the simulation if it crosses the obstacle bound- 
ary. In the analyzed run the simulation box contains 
about 12 million H + ions and about 0.2 million + ions 
at every time step. The average number of macroions per 
a cell is 30 and the timestep At = 0.04 s. 

In the model the interplanetary magnetic field is de- 
rived from the velocity field that represents the flow of 
a laminar incompressible fluid around a spherical obsta- 
cle by assuming that the magnetic field is "frozen in" to 
the laminar incompressible flow. The same approach has 
been used earlier in a 3-D QNH model [14]. If the flow at 
infinity is along the (negative) x-axis, the magnetic field 
outside of the obstacle can be shown to be 
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where Bq is the magnitude of the IMF and r = 
\/ x 2 + y 2 + z 2 . The radial component of the magnetic 
field on the surface of the obstacle is zero. The total 
magnetic field in the simulation box outside of the obsta- 
cle is therefore 



B(r,t) =B 



laminar 



(r)+$B(r,i), 



(9) 



where B^ aminar (r) is time independent magnetic field 
derived from Eq. |HI and 5B(r, t) is time depen- 
dent self-consistent magnetic field. It should be noted 




Figure 2: An example of magnetic field which is "frozen in" 
to the laminar incompressible flow around a sphere. Such a 
constant magnetic field is used in this work to generate the 
IMF x- and ^/-components. 



that B/ am ^ nar (r) is derived from the potential field 
(^laminar = V(/)) and therefore there are no electric 
currents associated with the B| aminar (r) field outside the 
obstacle (J ~ V x V(/) = 0). The magnetic field compo- 
nent ~Biaminar(Y) is associated with electric currents on 
the surface of the obstacle. 

Equations |H1 are needed to introduce the IMF x- 
component into the simulation box. In the previous 3-D 
QNH Venus model the IMF ^-component has been as- 
sumed to be zero 

urn 

which is not a good approxima- 
tion at Venus where the nominal Parker spiral angle is 
about 36° and, consequently, the x-component is often 
larger in magnitude that the y or z components. Imple- 
mentation of the IMF ^-component enables us to study 
its associated set of asymmetries. 

In this work Equations |H1 are used to model the full 
IMF, and not only its ^-component, by rotating the lam- 
inar flow around the z-axis in order to produce a mag- 
netic field which is equal to the IMF at infinity. Such 
a magnetic field is derived from Equations |H1 by rotating 
the magnetic field around the z-axis. Figure [21 shows an 
example of the magnetic fields used in this work. 



B. Plasma sources 

The simulations contain three H + sources: the solar 
wind, the cold atmospheric hydrogen and the hot neutral 
hydrogen exosphere. Likewise there are two + sources: 
the cold atmospheric oxygen and the hot neutral oxygen 
exosphere. The major source for H + is the solar wind. 
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In the simulation the solar wind H + ions are injected to 
the simulation box at the front face at x = 3Ry- The 
planetary hydrogen and oxygen densities are similar, or 
slightly simplified versions, of the density profiles that 
have been used recently to study generation of ENAs at 
Venus 

1. Hydrogen atoms 

The cold atmospheric hydrogen is modelled by the 
Chamberlain exosphere model 

n(r) = n e^ (1/ro - 1/r) C(^), (10) 

where no is the hydrogen density at the planetocentric 
distance r = r*o, C 18 Chamberlain's partition function 
and 

p = ^g^. (id 

Here G, My, fc# and T are the gravitational constant, 
planet's mass, Boltzmann's constant and neutral tem- 
perature, respectively. At noon we adopt the parameters 
n = 1.32 x 10 11 m~ 3 and T = 285 K and at the midnight 
n = 2.59 x 10 12 m~ 3 and T = 110 K. Both noon and 
midnight have ro = 170 km + Ry- The densities between 
noon and midnight are obtained by linear interpolation 
with respect to the solar zenith angle (SZA). Chamber- 
lain's partition function differs only slightly from unity 
within the simulation box and, therefore, it's value was 
taken to be unity. 

The hot hydrogen exosphere (corona) is modelled at 
noon, the terminator and midnight by a function 

n(r) = e air+a2+as/r , (12) 

where the coefficients en are given in Table |I] Linear 
interpolation in SZA is used to get densities between SZA 
= 0°, 90° and 180°. 



SZA (deg.) 


ai [km x ] 


a 2 


a 3 [km] 





-6.2625 x 10" 5 


15.4817 


3.6414 x 10 4 


90 


-8.4607 x 10" 5 


15.9944 


2.9743 x 10 4 


180 


-6.2309 x 10" 5 


15.2723 


4.3781 x 10 4 



Table I: The constants used in Equation El to model the hot 
hydrogen exosphere (corona). 



ions originating from the cold oxygen component is mod- 
elled by emitting oxygen ions from the obstacle bound- 
ary. The particle flux has its maximum at the noon and 
it decreases as cos(SZA) towards the terminator. At the 
nightside the oxygen ion emission rate is 10% of the noon 
value. The total + emission rate from the obstacle 
boundary is 10 25 s _1 . 

The hot oxygen exosphere is modelled by Equation ITU1 
with the following values: at noon no = 7.5 x 10 10 m -3 , 
r = 200 km and T = 6400 K and at midnight n 
2 x 10 9 m- 3 , r = 300 km and T = 4847 K. Values for 
other SZA values are obtained by linear interpolation. 

3. Ionization 

The hydrogen ion production rate g H + [s _1 m -3 ] and 
the oxygen ion production rate q Q + [s _1 m -3 ] in sunlight 
are modelled as 

«H+(r) = /h+W »H(r), (13) 
fo+(r) = /$(r)n (r), (14) 

where /^+(r)(/Q°+(r)) [s _1 ] is the total ionization fre- 
quency and nn(r) (no(r)) [m -3 ] the total density of the 
hydrogen (oxygen) atoms at r. and are zero 

within the optical shadow. The three most important 
ionization processes at Venus are photoionisation, elec- 
tron impact ionisation and charge exchange [see, for ex- 
ample, In this work only EUV ionization is taken 
into account for simplicity. Table ITT1 shows the used con- 
stant photoionisation frequencies for solar minimum con- 
ditions and the corresponding total ion EUV production 
rates within the simulation box. 



Source 


Ionization freq. 
/ (1/s) 


Total rate 

(#/s) 


Comments 


H hot — ► H + 

cold 

Ohot 0+ 

o cold -► 0+ 


1.39 x 10" 7 
4.55 x 10~ 7 


6.42 x 10 24 

4.09 x 10 24 
10 25 


EUV ionization 

EUV ionization 
From the obstacle 



Table II: Summary of H + and + production rates from the 
hot H neutral corona and cold hydrogen component and from 
the hot O neutral corona. The photoionization frequencies 
(/) are for the quiet solar EUV condition. Oxygen ions from 
the cold + component is modelled by emitting oxygen ions 
from the obstacle boundary. Note that part of the generated 
oxygen ions hit to the obstacle boundary and that the total 
escape rate from the simulation box is about 2.2 x 10 24 s _1 . 



2. Oxygen atoms 

The atmospheric scale height of the cold oxygen pop- 
ulation is much smaller than the grid size on the simu- 
lation (= 605 km) and therefore the emission of oxygen 



The presented run corresponds to the situation after 
240 seconds from the start of the simulation. At that 
time the simulation has reached a quasi-stationary state 
where the total energy within the simulation box does 
not increase much (Figure Et)- Also, at that time the 
total escape rate of + ions from the simulation box has 
reached a relatively stationary rate of about 2.2 x 10 24 s _1 
(Fig. Bb). 
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Tim a- jsl 

Figure 3: (a) The time evolution of the total energy density 
(kinetic energy density + magnetic energy density) inside the 
simulation box. (b) The number of escaping + ions from 
the side faces of the simulation box. 



III. RESULTS 

In this sections several plasma and magnetic field pa- 
rameters are presented on the XY and XZ planes. The 
used grid size of 0.1 Ry prevents us from a detailed study 
of the position and the shape of the various boundary 
layers. Nevertheless, a conic shape is included in the fig- 
ures to help the eye in catching asymmetries more easily 
and to simplify comparisons between various plasma and 
magnetic field parameters. The conical shape (l5l | 



r focus = (#0,0,0) (16) 

is used to provide a rough estimation for the shape of the 
bow shock. Here r is the radial distance from the focus 
point, r focus to the bow shock, e is the eccentricity, 9 is 
the angle between the a point of a bow shock and the pr- 
axis and L is the distance of the bow shock from the focus 
point on the x = xq plane. In this paper the conical shape 
on the XY and XZ planes is derived from Equations 
[El -[HD fc> r the parameters xo = O.QRy, e = 0.87 and 
L = l.SRy- 

Figure E| presents the density and bulk velocity of pro- 
tons in the XZ and XY planes. The proton density has 
increased and the velocity decreased at the bow shock. 
At the nightside behind the planet the density decreases 
again, becoming smaller than the undisturbed solar wind 
density. A clear asymmetry can be identified between 
the +Esw hemisphere (z > 0) and the -Esw hemisphere 
(z<0). 

It can be seen in Figs. Et andEfc that the bow shock 
is more blunt and further away from the planet on the 
+Esw hemisphere than on the -Esw hemisphere. Fur- 



thermore, the bulk velocity on the nightside in the tail 
is higher on the +Esw hemisphere than on the opposite 
hemisphere. Only a slight asymmetry can be found be- 
tween the BSperp hemisphere (y < 0) and on the opposite 
BS P ar hemisphere (y > 0). Note also that the distance 
of the bow shock from the planet is not identical in the 
XY and XZ planes. 

Figure [SI presents plasma parameters similar to those 
in Fig. E| but now for the oxygen ions. In this case a 
very clear +Esw/-Esw asymmetry is found. On the - 
Esw hemisphere a much sharper density gradient exists 
between the bow shock and the planet than on the +Esw 
hemisphere. Also, the bulk velocity of + ions near the 
bow shock and in the solar wind is much higher on the 
+Esw hemisphere than on the opposite hemisphere. In 
addition, the density of the oxygen ions near Venus is 
higher on the -Esw hemisphere than on the +Esw hemi- 
sphere. These asymmetries can be understood on the 
basis of the direction of the convective electric field and 
the large gyroradius of + ions compared to the size 
of the planet: the convective electric field accelerates 
newly born + ions away (against) from the planet on 
the +Esw(-Esw) hemisphere in the XZ plane. A slight 
BSperp/BSpar asymmetry arises in n(0 + ) the + density 
being slightly higher on the BS per p hemisphere than on 
the BSpar hemisphere. 

Figure El shows how the plasma parameters are asso- 
ciated with the properties of the magnetic field. The 
bow shock can be identified clearly from an enhanced to- 
tal magnetic field and enhanced B x component near the 
white line in Fig. and the black line in 6c, respectively. 
Increase of B x in the middle of the tail is associated with 
the magnetic tail lobe on the BS per p hemisphere, as will 
be seen later. 

In the XY plane a clear BS pe rp/BS par asymmetry can 
be found. The total magnetic field is much larger on 
the BSperp hemisphere (y < 0) than on the opposite 
hemisphere. In addition, the magnitude of B x on the 
magnetic lobe on the BS per p hemisphere is much larger 
than \B X \ on the opposite hemisphere. The asymmetry 
is associated with the IMF x-component: x-component 
of the ~Biaminar in Eq. EH is positive and thus it decreases 
\B X \ in the BS par hemisphere where B x is negative in the 
magnetic lobe. On the BS perp hemisphere B x is positive 
in the magnetic lobe and on that hemisphere ^laminar 
increases \B X \. 

Plasma and magnetic field in the Venusian tail is stud- 
ied in more detail in Fig. Note that the density of 
+ ions is in its maximum downstream of the bow shock 
in the same region where the proton density is low. The 
+Esw/-Esw asymmetry can be seen in all plasma param- 
eters: the bow shock is farther away from the planet in 
z > than in z < (Fig. OO, the proton velocity is 
higher in z > than in the opposite hemisphere (Fig. 
Cb), the maximum n(0 + ) is sited in the z < hemi- 
sphere and + ions can be found in the solar wind in the 
z > hemisphere (Fig. Et), the + ions obtain higher 
velocities in the +Esw hemisphere (z > 0) than in the 
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Figure 4: Properties of the H + ions near Venus: proton density [cm 3 ] on (a) the XZ plane and on (b) the XY plane and the 
bulk velocity [km s _1 ] on (c) the XZ plane and on (d) the XY plane. 




Figure 5: Properties of Venusian + ions: Particle density [log 10 ([cm 3 ])] (a) in XZ plane and (b) in XY plane and the bulk 
velocity [km s _1 ] (c) in XZ plane and (d) in XY plane. 



7 




Figure 6: Magnetic field near Venus: The total magnetic field [nT] (a) in XZ plane and (b) in XY plane and the x-component 
[nT] (c) in XZ plane and (d) in XY plane. 



opposite hemisphere (Fig. 01). The magnetic field is also 
asymmetric with respect to the direction of Esw- 

A somewhat unexpected feature in Figures and[Jf 
is the very clear asymmetry between the magnetic tail 
lobes. In the analyzed IMF direction case the electric 
currents in the tail results in the magnetic tail lobes hav- 
ing B x > at y < and B x < at y > 0, As already 
noted before, the IMF x-component, which is modelled 
by ^laminar in this study, decreases the total magnetic 
field at the y > hemisphere and increases the field at 
the opposite hemisphere. The magnetic tail lobe is there- 
fore much stronger at y < than at y > 0. 

A similar asymmetry can be found in the strength of 
the magnetic field between the y < and y > hemi- 
spheres. The magnetic field is clearly enhanced at the 
bow shock and in the magnetosheath in the BS perp hemi- 
sphere (y < 0). At the BS par hemisphere (y > 0), in- 
stead, the enhancement of the field is much weaker. In 
fact, practically no enhancement of the magnetic field 
can be seen in the BS par hemisphere at the z = plane. 

Figure |H1 illustrates how the morphology of the mag- 
netic field is associated with the BS perp /BS par asymme- 
try. On the BSperp hemisphere magnetic field piles up 
against the magnetic tail lobe, resulting in an enhanced 
magnetic field. On the BS perp hemisphere the sign of the 
B x at the magnetic tail lobe is the same as the sign of the 
B x in the solar wind. On the opposite BS par hemisphere 
the sign of the B x in the magnetic tail lobe is opposite 
to the sign of the B x in the solar wind and the direction 



of the field line therefore must change. 

IV. DISCUSSION 

This paper presents the first step to study the Venusian 
plasma environment with a new 3-D QNH model. To the 
authors' knowledge it is the first published QNH model 
study which uses a realistic Parker spiral IMF where the 
IMF x-component is larger in magnitude than the mag- 
netic field components perpendicular to the flow. 

The study suggests that the model can reproduce the 
basic observed plasma and magnetic field regions and 
boundaries near Venus, implying a future potential of 
the developed numerical approach. However, the study 
leaves many open issues for future studies. The grid res- 
olution is not fine enough to resolve the position and the 
shape of the bow shock and the inner structure of the 
magnetic barrier. Especially, the grid resolution makes 
it impossible to include a self-consistent ionosphere. Fur- 
thermore, at the present stage of the model the total + 
emission rate from the obstacle boundary is a freely cho- 
sen parameter. In the future, its value can be adjusted by 
comparing the simulated + densities with the observed 
densities. Similar approach has been used earlier to esti- 
mate the total ion escape rate at Mars and a similar type 
of estimation can be made when direct + measurements 
from ASPERA-4/ Venus Express mission becomes avail- 
able in 2006. In the future one of the biggest challenges 
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Figure 7: Plasma and the magnetic field at x = —2Ry plane: 
in km/s -1 ] of protons, (c) the density [n(0 + ) in log 10 ([cm -3 ]) 
planetary + ions, (e) the total magnetic field [\B\ in nT], and 



in modelling is to implement a more realistic ionosphere 
to the 3-D QNH model. 

The second limitation of the presented study is that is 
presents the solution only for one set of upstream param- 
eters. In practice, Venus is under the influence of varying 
plasma and magnetic field conditions [see, for example, 
How the model responds to the upstream parame- 
ters, which also have time dependence, will be a topic of 
future studies. 

It is worth recalling that in the presented run the re- 
sistivity is constant outside Venus. The used resistivity 
was not 'optimized' to its minimum value and the value 
is artificial in the sense that its value can be smaller in re- 
ality. The specific resistivity value was adopted because 
it has been found to be large enough to be used for var- 
ious upstream parameter runs (not shown in this paper) 
to reduce fluctuations formed near the obstacle bound- 
ary (Venus). The role of the non-uniform, for example, 
spherically symmetric resistivity model, will be studied 
in the future when a smaller grid size near Venus than 
used in this paper will be used. 

It is also worth noting that there is a large mag- 
netic field data set available from Pioneer Venus Orbiter 




(a) the density [n(H + ) in cm -3 ] and (b) the velocity [v(H + ) 
scale] and (d) the velocity [v(0 + ) in km s _1 ] of the escaping 
(f) the x-component of the magnetic field [B x in nT]. 



(PVO) over ten year measurements. The data set makes 
it possible to make a detailed quantitative comparison be- 
tween the developed QNH model and the observations. 
In fact, preliminary comparisons between the model and 
the magnetic field measurements made by PVO when the 
spacecraft clearly crossed both magnetic tail lobes [seeHH 
have shown that the developed Venus model is capable 
of reproducing many of the observed features (figures not 
shown). In the model the high velocity + ions were also 
found on the +Esw hemisphere (c.f. Fig. Efc), that is, on 
the same hemisphere where PVO observed fast moving 
escaping + ions, for example, in the far tail [see, for 
example, [Hj. A detailed comparison between PVO mag- 
netic field and ion observations and the model will be a 
topic of future studies. 

It is interesting to note that the BS pe rp/BS par asym- 
metry of the bow shock on the XY plane (Fig. EJd and 
@Ji) is in agreement with observations [seeQi|- One may 
anticipate that a detailed comparison between the data 
and the model will provide a new insight about the ap- 
plicability of the global 3-D QNH model to reproduce 
plasma parameters at various upstream conditions. The 
developed QNH model can be used to separate effects of 
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Figure 8: Draped magnetic field lines at Venus. Note the 
asymmetry between the quasi-perpendicular bow shock hemi- 
sphere, BSperp (y < 0), and the quasi-parallel bow shock 
hemisphere, BS pa r (y > 0). 



different kinds. In some sense the model can also be re- 
garded to provide a possibility to "filter" away from the 
measurements such effects that result from the very basic 
laws of nature (Lorentz force, Faraday's and Ampere's 



law and the conservation of the electron's momentum) 
under the used boundary conditions from effects that 
possibly result from unknown factors not implemented 
in the model. For example, the role of the electron im- 
pact ionization and the role of the ion-neutral collisions 
remains to be studied, especially, which kind of global 
effects the omitted process may cause. 

Finally, possible artefacts result from the adopted mag- 
netic field model (Equations |B| have to be studied by 
comparing runs made for different magnetic field config- 
urations. 



V. SUMMARY 

A three dimensional self-consistent numerical model 
is developed to study the Venus-solar wind interaction. 
An analysis of the runs made for a Parker spiral IMF 
direction suggests that several asymmetries take place 
at Venus. The observed +Esw/-Esw asymmetry re- 
sults from ion finite gyroradius effects and from the 
escaping + ions while the IMF ^-components caused 
BS perp /BS P ar asymmetry. The density of the escaping 
planetary + ions in the magnetosheath and in the so- 
lar wind is found to be higher in the hemisphere where 
the convective electric field points in the undisturbed so- 
lar wind than in the opposite hemisphere. The analysis 
implies that the properties of plasma and magnetic field 
at given points near Venus can be anticipated to vary 
noticeably in concert with varying upstream parameters. 
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